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Abstract 

We analyze the collective behavior of neutrinos and antineutrinos in a dense back- 
ground. Using the Wigner transform technique, it is shown that the interaction can 
be modelled by a coupled system of nonlinear Vlasov-like equations. From these 
equations, we derive a dispersion relation for neutrino-antineutrino interactions on 
a general background. The dispersion relation admits a novel modulational instabil- 
ity. The results are examined, together with a numerical example, and we discuss the 
induced density inhomogeneities using parameters relevant to the early Universe. 

PACS numbers: 13.15.+g, 14.60Lm, 97.10Cv, 97.60Bw 



1 Introduction 



Neutrinos have fascinated people ever since they were first introduced by Wolf- 
gang Pauli in 1931. Since then, neutrinos have gone from hypothetical to an 
extremely promising tool for analysing astrophysical events, and neutrino cos- 
mology is now one of the hottest topics in modern time due to the discovery 
that neutrinos may be massive [1]. Because of its weak interaction with other 
particles, neutrinos can travel great distances without being affected appre- 
ciably by material obstacles. They can therefore give us detailed information 
about events taking place deep within, e.g. supernovae. Furthermore, since 
the neutrinos decoupled from matter at a redshift z of the order 10 10 , as 
compared to z ~ 10 3 for photons, it is possible that neutrinos could give 
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us a detailed understanding of the early Universe, if such a signal could be 
detected [2]. Massive neutrinos have also been a possible candidate for hot 
dark matter necessary for explaining certain cosmological observations, such 
as rotation curves of spiral galaxies [3]. Thus, massive neutrinos could have 
a profound influence on the evolution of our Universe. Unfortunately, due to 
the Tremaine-Gunn bound [4], the necessary mass of the missing particles (if 
fermions) for explaining the formation of dwarf galaxies seems to make neu- 
trinos of any species unlikely single candidates for dark matter. As a remedy 
to this problem, interacting hot dark matter has been suggested [5,6], since 
the interaction prevents the free-streaming smoothing of small scale neutrino 
inhomogeneities. Thus, dark matter in astrophysics is not only a mystery but 
it also plays an essential role in determining the dynamics of the universe, 
its large scale structures, the galaxies and superclusters. However, so far, the 
suggested "sticky" neutrino models have not been successful in dealing with 
the dwarf galaxy problem [5]. 

A first successful indication that neutrinos have a non-zero mass came in 1998 
through laboratory experiments of atmospheric neutrinos and their oscillations 
[7]. Although the allowed neutrino masses encompass a wide range 1 , it is 
currently believed that neutrinos have masses below 2 eV. This conclusion 
is furthermore supported by independent cosmological observations (see, e.g., 
[9]). Thus, the masses of neutrinos are indeed very small, and the classical 
analysis by Tremaine and Gunn would thereby indicate that neutrinos can 
in no way be considered as a sole candidate for dark matter. This conclusion 
will in this paper be re-analyzed within the electro-weak framework, where 
neutrino-neutrino interactions occur as a natural consequence of the theory. 

Thus, in this Letter, we consider the nonlinear interaction between neutri- 
nos and antineutrinos in the lepton plasma of the early Universe, adopting 
a semi-classical model. Neutrinos and antineutrinos interact with dense plas- 
mas through the charged and neutral weak currents arising from the Fermi 
weak nuclear interaction forces. Charged weak currents involve the exchange 
of the charged vector bosons associated with the processes involving interac- 
tions between leptons and neutrinos of the same flavor, while neutrino weak 
currents involve the exchange of the neutral vector bosons associated with 
processes involving neutrinos of all types interacting with arbitrary charged 
and neutral particles. Asymmetric flows of neutrinos and antineutrinos in the 
early Universe plasma may be created by the ponderomotive force of nonuni- 
form intense photon beams or by shock waves. Here, using an effective field 
theory approach, a system of coupled Wigner-Moyal equations for nonlinearly 
interacting neutrinos and antineutrinos is derived, and it is shown that these 
equations admit a modulational instability. Finally, we discuss the relevance 



1 Some estimates even support the notion that neutrinos may contribute up to 20% 
of the matter density of the Universe [8] . 
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of our results in the context of the dark matter problem, and it is moreover 
suggested that the nonlinearly excited fluctuations could be used as a starting 
point for obtaining a better understanding of the process of galaxy formation. 
It turns out that the short-time evolution of the primordial neutrino plasma 
medium in the temperature range 1 MeV < T < 10 MeV is governed by colli- 
sionless collective effects involving relativistic neutrinos and antineutrinos. 



2 Dispersion relation and the motion of neutrino bunches 

As a primer, we will study the implication of the known dispersion of neutrinos 
on a thermal neutrino/anti- neutrino background, using the eikonal represen- 
tation and the WKBJ approximation. 

Suppose that a single neutrino (or anti-neutrino) moves in a Fermionic sea 
composed of neutrino-antineutrino admixture. The energy E of the neutrino 
(antineutrino) is then given by (see, e.g. [10,11]) 

E = ^fp 2 c 2 + m 2 c 4 + V±(r, t), (1) 

where p is the neutrino (anti-neutrino) momentum, c the speed of light in vac- 
uum, and m the neutrino mass. The effective potential for a neutrino moving 
on a background of it's own flavor and in thermal equilibrium is given by 2 
[10] (see also [12,13,14,15]) 

n — n), (2a) 

while the potential for a neutrino moving on a background of a different flavor 
is 

V±(r, t) = ±V2G F (n - n), (2b) 

where G F /(hc) 3 rs 1.2 x 10~ 5 GeV" 2 , G F is the Fermi constant, n (n) is the 
density of the background neutrinos (antineutrinos), and + (— ) represents 
the propagating neutrino (antineutrino). Expressions (2) are valid in the rest 
frame of the background. As seen from (1) and (2), while neutrinos moving 
in a background of neutrinos and antineutrinos change their energy by an 
amount ~ G F (n — n), the antineutrinos change their energy by ~ — G F {n — n) 
[16]. The extra factor of 2 in (2a) as compared to expression (2b) comes from 
exchange effects between identical particles [13]. 

The relation (1) can be interpreted as a dispersion relation for relativistic and 



2 For a more detailed description of the potential, see the next section. 
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nonrelativistic neutrinos, with the identifications E = Tiu and p = hk, i. e. 



where % is the Planck constant divided by 27r. By using the eikonal represen- 
tation (viz. E — > Tiuoq — ihd/dt and p — > hko + ihV) and the WKBJ approx- 
imation [17,18] (viz. \d&/dt\ < uol&l and |W| < l^oll^l), we obtain from 
Eq. (3) a Schrodinger equation for slowly varying modulated (by long-scale 
density fluctuations) neutrino (anti-neutrino) wave function &(r,t) (i.e., neu- 
trino bunches) (see also Ref. [16] for a similar treatment of neutrino-electron 
interactions) 

<P - = 0, (4) 

where v g = ck (k 2 + m 2 c 2 /h 2 )" 1 ^ 2 is the group velocity 3 of relativistic neutri- 
nos and antineutrinos which have similar energy spectra, 7 = (1 — t^/c 2 ) -1 / 2 
is the relativistic gamma factor, n = k /\k \, and k is the vacuum wavevec- 
tor. Suppose now that the neutrino bunches themselves are nearly in thermal 
equilibrium (to be quantified in the next section). Then, we have the case of 
self-interacting neutrinos and anti-neutrinos, and the densities in the potential 
V± is given in terms of the sums 



d_ 
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+ « • v hp + 
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i=i 



i=l 



i=l 



i=l 



where <Pi + and !Pj_ are the neutrino and antineutrino wave functions (with % 
numbering the wave functions), respectively, and the angular bracket denotes 
the ensemble average. In this case, the relativistic neutrino and antineutrino 
wave packets are comoving with the background, and Eq. (4) thus yields 

where Vj_ = V 2 - (n • V) 2 , V 2 = (n • V) 2 . Expressions (2a) and (5) 
reveal that self-interactions between relativistic neutrinos and antineutrinos 
produce a nonlinear asymmetric potential in Eq. (6). By further rescaling the 



3 We note that when the scalelength of the density inhomogeneity is comparable 
to the wavelength of the modulated neutrino wave packets, we must modify the 
coupled Schrodinger equations to account for differing group velocities of neutrinos 
and antineutrinos in a Fermionic sea. We expect a shift in the momentum of Eq. (13) 
and a slower growth rate of the modulational instability of neutrino quasi-particles 
involving short-scale density inhomogeneities. 
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coordinate along n , Eq. (4) can finally be written as the coupled system 



i-gp + \ V 2 ^ ± =F P(n - n)!Pi± = 0, (7) 

where a = h/rwy, and /3 = 2v / 2G\f/7i for neutrinos moving on the same flavor 
background. 

Equation (7) shows that this approach can lead to some interesting effects. 
The case of a single self-interacting neutrino bunch shows that the formation of 
dark solitary structures is possible. Furthermore, the slightly more complicated 
case of two interacting bunches, either of the neutrino-neutrino or neutrino- 
anti-neutrino type, can result in splitting and focusing of the wave packets 
[19]. 



3 Kinetic description 



In the preceding section, we investigated the case of a neutrino bunch close 
to thermal equilibrium. In general, this may of course not be the case, and 
Eq. (2) must be modified. The more precise form of the potential V± for equal 
species due to neutrino forward scattering is given by [20] 



V±(t, r,p, f i± ) = ±2V2G F Jdq(l-p-q) 



■ M N 

12 /*+(*> - 12 /*-(*> r > 9) 

i=l i=l 



(8) 

where hatted quantities denote the corresponding unit vector, and fi+(t, r, q) 
(fi-(t,r,q)) is the neutrino (anti-neutrino) distribution function correspond- 
ing to bunch %. The distribution functions are defined to be normalized such 
that 

Ui(t,r) = f dqf i+ (t,r,q), ni(t,r) = J dq fi-(t,r, q), (9) 
where Hi (n«) is the number density of the i th neutrino (anti-neutrino) bunch. 

The first thing to notice is that when the distribution is thermal, the po- 
tential (8) reduces exactly to (2a). Secondly, when the neutrinos have an al- 
most thermal distribution, i.e. the corresponding distribution function may 
be expressed as (dropping the indices for notational simplicity) f(t,r,p) = 
fo(p) + Sf(t,r,p), where \Sf\ |/ |, we obtain the following form of the 
potential 

r , / m n 

V ± (t,r,p;f l± ) = ±2V2G F 



(n-n)- J dq(p- q) I Sf i+ - ^ 5f\ 



i=i i=i 



(10) 
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The last term is small and may therefore be neglected, and we obtain V±(t, r) « 
accordance with expressions (2a), thus justifying the equa- 
tion of motion (7). 

Now, we define a distribution function for the neutrino states by Fourier trans- 
forming the two-point correlation function for <?±, according to [21] 

fi±(t,r,p) = J dye^y/ n (*: ± (t,r + y/2)V l± (t,r-y/2)), (11) 

where p represents the momentum of the neutrino (antineutrino) quasi-particles 
(note that the ensemble average was not present in the original definition [21], 
but has important consequences when the phase of the wave function has a 
random component [22]; for similar treatments of optical beams and quantum 
plasmas, see [22] and [23], respectively). We note that with the definition (11), 
the following relation holds 

m±(t,r)\ 2 } = j dpf i± (t,r,p). (12) 

Thus, using (11) and (6) together with the potential (8) we obtain the gener- 
alized Wigner-Moyal equation 

dfi± | p df i± 2\4 
dt m7 dr h 

for fi±, where the sin-operator is defined in terms of its Taylor expansion, 
and the arrows denote the direction of operation. In the case of the potential 
(2a), the last term in the sin-operator drops out, and Eq. (13) reduces to the 
standard Wigner-Moyal equation [21]. 

Retaining only the lowest order terms in h (i.e. taking the long wavelength 
limit), we obtain the coupled Vlasov equations 

' d_ (p_ dV±\ d_ 
dt dp J dr 

The term dV±/dp represents the group velocity. While higher order group 
velocity dispersion is present in (13), this is not the case in (14). Thus, in- 
formation is partially lost by using the Eq. (14). Furthermore, while Eq. (14) 
preserves the number of quasi-particles, Eq. (13) shows that this conclusion 
is in general not true, i.e. the particle number in a phase-space volume is not 
constant, and the higher order terms d n V±/dr n may moreover contain vital 
short wavelength information. Equations similar to (14) have been used to 
study neutrino-electron interactions in astrophysical contexts [11]. 

Suppose now that we have small amplitude perturbations on a background of 



hid d d d 
2 \ dr dp dp dr 



fi± = 0, (13) 



dV± dfe 
dr dp 



= 0. 



(14) 
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constant neutrino and antineutrino densities rii = riio and rii = n^o, respec- 
tively, i.e. 

fi±(t, r, p) = f i0 ±(p) + 8f i± (p) exp[t(K ■ r - f2t)], (15) 

and \5fi±\ <C |/io±|, where K and i? is the perturbation wavevector and fre- 
quency, respectively. Thus, Eqs. (13) give 



p K 

m7 




i± 



2 *r ih „ d 

+ ^ 111 \ ~2 dp 



f, 



i0± 



0, 



(16) 



where 5V± = V±(t,r,p; Sfi±) and V ± = V±(t,r,p; f i0 ±) from Eq. (8). Elimi- 
nating 5fi± from (16), using SV- = — we have 



SV+(P) = J d q{l-p.q) SV + (q) 



x 



M 

E 



sin (f K ■ V) Mq) 



rn-y 



N 

+E 



sin 



(17) 



Assuming that Sfi±(p) is a symmetric function of p implies that SV± is inde- 
pendent of p, and Eq. (17) simplifies to the dispersion relation 



1 = 



Ay/2 id 



N 



+E 



M 

E 



sin(fK.|)/, 0+ (q) 



j^-^-fsKf^-£)H + (<?) 



sin 



K 



d_ 
dq 



fio-(q) 



(18) 



where we have dropped the arrows indicating the direction of operation. Note 
that if the background distribution is thermal, V ± is independent of p, and 
the last term in the denominators of Eq. (18) vanishes. 
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3.1 The one- dimensional case 



The simplest way to analyse the dispersion relation (18) is to reduce the dimen- 
sionality of the problem. We therefore first look at the one- dimensional case, 
where we may use the identity 2 sin (4r J;) h(p) = i \h(p + ^) - h(p - ¥f) , 
in order to rewrite the dispersion relation (18) as 




hK/2)- f l0+ (q- hK/2) 
Q - qK/rwy + A+(q) 



v /,o_ (g + UK 1 2) - f^jq - hK/2) \ 
+ Q-qK/m^ + A^q) J ' 1 ' 

where we have introduced A±(q) = [Vo±(? + hK/2) — V ±(q — kK/2]/h. In 
the case of mono-energetic beams, i.e. f i0+ (p) = nio5(p — p i0 ) and fio~(p) — 
nioS(p —pio), Eq. (19) reduces to 



1 = 



2V2Gj 



M 
i=l 



m7 



+ A + (p i0 + hK/2) - A + (p i0 - hK/2) 



x 
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'hK 2 \ 2 
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K ) [A+feo + + A + ( Pl0 - hK/2)} 
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hK 2 

2my 



[A + (p l0 + hK/2) - A + ( Pl0 - hK/2)] 
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+A+(p l0 + hK/2)A + (p l0 - hK/2) 
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+ A_(p i0 + hK/2) - A_(p i0 - hK/2) 



hK 2 



2mj 



[A_fe + hK/2) + A_fe - hK/2)] 



hK 2 
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[A_(p l0 + hK/2) - A_fe - hK/2)} 



+A_(p l0 + hK/2)A_(p i0 - hK/2) 



(20) 
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where 

V ±(p) = ±2V2G F 
by Eq. (8). 



M 



N 



(n - no) - sgn(p) ^ n i0 sgn(p i0 ) - n *o sgn(p i0 ) 



Vi=l 



i=l 



(21) 



Let us look at the simplest case of interacting neutrinos and antineutrinos 
with M — N — 1. We assume that they have equal densities no = no, and are 
counter-propagating, i.e. p = — Po > 0. From (21) we obtain the potential 



Vb±(p) = =F4v^Gj?sgn(p)7io 



while Eq. (20) yields 
(UK 2 



+ 



\ m7 



— 2ae < 



m7 , 



ft - 
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(22) 



p ^\ ^ 



1 1 m7 J 2m7 J 



m7 y 2m7 ; 



1, (23) 



where a = 2\ / 2GFn /h and e = 1 — sgn(p — ft-K") = 0, 1, or 2 when p > ft-K", 
Po = ftK, or < ^-K", respectively. Thus, for the case e — 0, the growth rate 
is given by (see Figs. 1 and 2) 



£!. 



Av 2 



2m7 



+ Av 2 v 2 F + Vp — v 2 — Vp — 



'HIT 
2m7 



(24) 



where r = —%Q is the instability growth rate, and vp = 2\/2 GVno/nry. Note 
also that, as expected, in the limit v — > 0, the instability disappears, just 
stating the well-known fact that there must be a non-zero relative velocity 
between the beams in order for the instability to occur. As F 2 is positive, we 
have 

"2< (tztt— ) < (-) (25) 



Vp 



2m r )Vp 



Vp 



i.e. 



l <£ <£ (l-2v 2 p/v 2 )- 1/2 , (26) 

where we have introduced the length scales £ = 2tt/K and £ = h/{2rrvyv). 
Thus, a higher neutrino momentum can retain a smaller instability length 
scale. It is clear from (24) that (i) the instability will remain for arbitrary 
velocities (see Figs. 1 and 2), and that (ii) the higher the neutrino velocity, 
the smaller the corresponding instability length scale £. 



-l 



+ 
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0.5 1 1.5 , 2 2.5 3 



Fig. 1. A contour plot of the values of (r/Kv) 2 , when e = 0, for which the instability 
occurs. The function (r/Kv) 2 is constant along the contours, and plotted in terms 
of the variables (vf/v) 2 and (h/2vm^i) 2 . Outside the contours F 2 < 0. 




Fig. 2. The same plot as in Fig. 1, but for e = 2. 
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3.2 Partial incoherence and thermal effects 



Partial incoherence can in general lead to lower growth rate, similar to inverse 
Landau damping. As an example of the results of stochastic effects, e.g. ther- 
mal fluctuations, we look at the following example. Let the indeterminacy of 
the wave function manifest itself in a random phase tp(x) of the background 
wave function, with the width Ap defined according to {^~ l< p{ x +y/ 2 ) q^^-v/ 2 )^ — 
e -Ap|j/|/ft_ rj ue to thjg ranc [ om spread, the modulational instability will be 
damped, as will be shown below. The Wigner function corresponding to the 
random phase assumption is given by the Lorentz distribution 

71" VP -PO) + AP 

With this, we obtain Eq. (24) with r — > _Td + ApK/rwy, where J 1 /) is the 
reduced growth rate. Thus, we see that the broadening tends to suppress the 
growth. Moreover, a positive growth rate r D requires 

2VAP < 1 

"XT < V (28) 

where r is given by Eq. (24). Hence, the general property of a spread in 
momentum space, here exemplified by a random phase, is to put bounds on 
the modulational instability length scale t. 

Incoherent effects among the neutrinos and anti-neutrinos can also be ap- 
proached for a background obeying Fermi-Dirac statistics, i.e. 



fo±{P) 



\n(4)k B T z 



± 



—GST 



(29) 



where we set M — N — 1, and assume n = n . Here, we have neglected the 
mass of the neutrinos (which will give us the correct result to lowest order). 
We will for simplicity assume that T± = T, so that the dispersion relation (19) 
takes the form 



1 = - ^) C ^T° II dp {n " P R / m ^ 1 [(! + ex P( c ^ + KKW/ksT))- 1 
+ (l + exp(c\p-hK/2\/k B T)y 1 ]. (30) 

The dispersion relation (30) cannot be solved analytically, but it can be ex- 
pressed according to 



1 = -Q [P(I(f2 n , K n )) + mg{Q n , K n )\ , (31) 
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where P(I(f2 n , K n )) is the principal value of the integral 

Q + K 2 Q - K 2 



1= / ^(l + e*)" 1 
Jo 

and g = g+ + g~, where 



+ 



(f2 n + K 2 ) 2 - K 2 x 2 ((2 n - K 2 ) 2 - K 2 n x 



(32) 



Q ±K 2 

g±(f2 n , K n ) = " " = , (33) 

Q = (A/\nA)(2^/2G F n /k B T)(m-ic 2 /k B T), and we have introduced the di- 
mensionless variables Q n = (hrwyc 2 / \k B T) 2 )fl and fT n = (hc/y/2 k B T)K. 
The constant Q gives the ratio of the potential energy contribution of the 
background and the individual neutrino energy to the thermal energy of the 
background. Furthermore, m r )c 2 ~ k B T, thus simplifying the expression for 
Q. The contributions from real and imaginary parts to the dispersion rela- 
tion are plotted in Figs. 3 and 4, respectively. Note that for very short length 
scales, i.e. large K, the quantity f2 n — K 2 becomes negative, and the imaginary 
part in Eq. (31) change sign, something which will not show using the long 
wavelength limit equation (14). This "quantum" behaviour can in principle 
lead to growth instead of damping of the perturbations (see Ref. [25] for a 
general discussion of this behaviour). We can obtain a quantitative measure 
of the growth rate as follows. For any fixed K n o, the dimensionless growth rate 
r n = —i Im Q n may be expressed as, denoting the value at i?„o by 0, 

r _ ^ (Q- 1 + PCjo)) (dg/dOno) - g (dP(I)/dn n0 ) 
n 71 n 2 (dg/df2 n0 ) 2 + (dP(I)/dn n0 ) 2 

to first order around (f2 n o, K n0 ). Thus, F n > if (Q _1 + P(iio)) (dlng/df2 n0 ) > 
dP{I)/dQ n Q. Moreover, using values given in Sec. 4, one can show that Q^ 1 ~ 
3x 10 9 . Thus, over a wide range of (i? n , K n ), Q^ 1 dominates the contribution to 
the growth rate, and a positive growth rate is implied as long as dg/df2 n0 > 0. 



4 Applications 



As a model for hot dark matter, massive neutrinos have for some time been 
one of the prime candidates, but as such they have faced the problem of the 
scale of the inhomogeneities they can support. Due to the conservation of 
phase-space density, the Tremaine-Gunn limit constrains the neutrino mass 
for isothermal spheres of a given size. For dwarf galaxies, for which there 
are ample evidence of dark matter [26], the necessary mass of the neutrino 
is uncomfortably large [4,27]. On the other hand, as pointed out by Raffelt 
& Silk [5], interacting dark matter can in principle change this picture. Here 
we see from Eq. (26) that as the neutrino momentum increases, the typical 
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2 4 6 8 10 

Fig. 3. A contour plot of the Cauchy principal value P(/(J7 n , K n )) as a function of the 
dimensionless variables J? n and K n . We note that P(I(f2 n ,K n )) > 0, being largest 
for small Q n and K n , and approaching zero at infinity. The uppermost contour has 

p(i(n n ,K n )) = o. 




2 4 6 8 10 



Fig. 4. A contour plot of the contribution g(f2 n , K n ) of the poles to the integral (30) 
as a function of the dimensionless variables Q n and K n . The darker areas represent 
negative values, the lighter positive values, and g is zero on the contour emanating 
from {a n ,K n ) = (0,1.8). 

length scale I of the inhomogeneity that can be supported by the modulational 
instability decreases. From the definition of £ , we note that as v tends to c, 
£ — > 0, and due to Eq. (26) the allowed scale of inhomogeneity becomes 
squeezed between two small values. On the other hand, if v ~ vp (a condition 
stating that the neutrino number density must reach extreme values), the 
upper inhomogeneity scale limit diverges. A minimum requirement for the 
effect to be of importance is that the instability growth rate is larger than 
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the Hubble parameter H . An estimate of the growth rate can be obtained as 
follows. At the onset of "free streaming" of neutrinos (i.e. their decoupling 
from matter and radiation) at z ~ 10 10 , the neutrino number density can be 
estimated as n ~ 2.1 x 10 38 m~ 3 (see, e.g. [28]). Furthermore, we assume that 
the neutrino mass is in the range m ~ 1 eV, and find vf — 2.1 x lO 6 ^ 1 / 2 m/s. 
The temperature of the neutrinos, given by T„ = (4/ll) 1//3 To(l + z) (To being 
the present day CMB temperature) [28], at neutrino decoupling is T„ ~ 2 x 
10 10 K. Thus, the thermal energy is roughly five orders of magnitude greater 
than the assumed rest mass of the neutrino, and in this sense the neutrinos 
can be well approximated as ultra-relativistic. In this case, using values of 
(UK / 2m r yv f) 2 in the middle range of the inequality (25), Eq. (24) gives r ~ 
2^/2GpTio/h ~ 16 x 10 10 s _1 for the values specified above. Assuming a critical 
density for the Universe, the Hubble time becomes H^ 1 ~ Hq 1 ^ + z)~ 3 / 2 ~ 
5 x 10 2 s at a redshift 10 10 , and thus r/H > 1. 

Although the two-stream instability may seem contrived as a cosmological ap- 
plication, the important issue displayed by this example is the non-gravitational 
growth of inhomogeneities, given a small perturbation of a homogeneous, al- 
though anisotropic, background. The fact that the growth rate exceeds the 
inverse of the Hubble time by many orders of magnitude makes it clear that 
the mechanism may be of some importance. Moreover, the analogous estimate 
for the Fermi-Dirac background, although done in a simplistic manner, in- 
dicates that the growth of the large K perturbations may be of importance. 
Note that this effect is a result of the use of the full Wigner-Moyal system, 
as compared to the Vlasov system (14), where these short wavelength effects 
are manifestly neglected. Furthermore, it could also be of interest to use the 
current formalism as a tool to investigate neutrino interactions within super- 
novae, where the two-stream instability scenario may occur as a more natural 
ingredient than perhaps within cosmology. 



5 Conclusion 

In conclusion, we have considered the nonlinear coupling between neutrinos 
and anti-neutrinos in a dense plasma. It is found that their interactions are 
governed by a system of Wigner-Moyal equations, which admit a modulational 
instability of the neutrino/antineutrino beams against large scale (in compar- 
ison with the neutrino wavelength) density fluctuations. Physically, instability 
arises because interpenetrating neutrino and antineutrino beams are like quasi- 
particles, carrying free energy which can be coupled to inhomogeneities due 
to a resonant quasiparticle-wave interaction that is similar to a Cherenkov 
interaction. Nonlinearly excited density fluctuations can be associated with 
the background inhomogeneity of the early Universe, and possibly counteract 
the free streaming smoothing of the small scale primordial fluctuations, thus 
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making massive neutrinos plausible as a candidate for hot dark matter. 
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